Astronomy & Astrophysics manuscript no. 5173 


© ESO 2008 


February 5, 2008 





Swift observations of GRB 050904: the most distant cosmic 

explosion ever observed. 

so ■ 

q ■ G. Cusumano 1 , V. Mangano 1 , G. Chincarini 2 3 , A. Panaitescu 4 , D.N. Burrows 5 , V. La Parola 1 , T. Sakamoto 6 ' 7 , S. 
Campana 2 , T. Mineo 1 , G. Tagliaferri 2 , L. Angelini 6 , S.D. Barthelmy 6 , A.P. Beardmore 8 , P.T. Boyd 6 , L.R. 
Cominsky 9 , C. Gronwall 5 , E.E. Fenimore 4 , N. Gehrels 6 , P. Giommi 10 , M. Goad 8 , K. Hurley 11 , S. Immler 6 , J.A. 
q \ Kennea 5 , K.O. Mason 12 , F. Marshal 6 , P. Meszaros 5 - 13 , J.A. Nousek 5 , J.P Osborne 8 , D.M. Palmer 4 , P.W.A. Roming 5 , 

A. Wells 8 , N.E. White 6 , B. Zhang 14 



\ 1 INAF-IASFPA, Via Ugo La Malfa 153, 90146 Palermo, Italy 
2 INAF - Osservatorio Astronomico di Brera, Via Bianchi 46, 23807 Merate Italy 
^ ' 3 Universita degli studi di Milano-Bicocca, Dip. di Fisica, Piazza delle Scienze 3, 1-20126 Milan, Italy, 
O ! 4 Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA 

C ■ 5 Department of Astronomy & Astrophysics, Pennsylvania State University,University Park, PA 16802 USA 
]Q • 6 NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 

National Research Council, 2001 Constitution Avenue, NW, TJ21 14, Washington, DC 20418, USA 
Department of Physics and Astronomy, University of Leicester, University Road, Leicester LEI 7RH, UK 
Department of Physics and Astronomy, Sonoma State University, Rohnert Park, CA 94928-3609, USA 
ASI Science Data Center, via Galileo Galilei, 00044 Frascati, Italy 
UC Berkeley Space Sciences Laboratory, Berkeley, CA 94720-7450 
Q-f 12 MSSL, University College London, Holmbury St. Mary, Dorking, RH5 6NT Surrey, UK 
q . 13 Department of Physics, Pennsylvania State University, PA 16802, USA 

$-1 ' 14 Department of Physics, University of Nevada, Box 454002, Las Vegas, NV 89154-4002, USA 

<Z3 ■ 

Received: ; accepted: 



o 

SO 

o 



7 
8 
9 
10 
1 1 



X 



ABSTRACT 



Context. Swift discovered the high redshift (z=6.29) GRB 050904 with the Burst Alert Telescope (BAT) and began observing with its narrow 
field instruments 161 s after the burst onset. This gamma-ray burst is the most distant cosmic explosion ever observed. Because of its high 
redshift, the X-ray Telescope (XRT) and BAT simultaneous observations provide 4 orders of magnitude of spectral coverage (0.2-150 keV; 
1.4-1090 keV in the source rest frame) at a very early source-frame time (22 s). The X-ray emission was monitored by the XRT up to 10 days 
after the burst. 

Aims. We present the analysis of BAT and XRT observations of GRB 050904 and a complete description of its high energy phenomenology. 
Methods. We performed time resolved spectral analysis and light curve modeling. 

Results. GRB 050904 was a long, multi-peaked, bright GRB with strong variability during its entire evolution. The light curve observed by 
the XRT is characterized by the presence of a long flaring activity lasting up to 1-2 hours after the burst onset in the burst rest frame, with 
no evidence of a smooth power-law decay following the prompt emission as seen in other GRBs. However, the BAT tail extrapolated to the 
XRT band joins the XRT early light curve and the overall behavior resembles that of a very long GRB prompt. The spectral energy distribution 
softens with time, with the photon index decreasing from —1.2 during the BAT observation to -1.9 at the end of the XRT observation. The 
dips of the late X-ray flares may be consistent with an underlying X-ray emission arising from the forward shock and with the properties of the 
optical afterglow reported by Tagliaferri et al. (2005b). 

Conclusions. We interpret the BAT and XRT data as a single continuous observation of the prompt emission from a very long GRB. The 
peculiarities observed in GRB 050904 could be due to its origin within one of the first star-forming regions in the Universe; very low metallicities 
of the progenitor at these epochs may provide an explanation. 

Key words, gamma rays: bursts, gamma-ray bursts: individual (GRB 050904) 



1. Introduction 



Send offprint requests to: G. Cusumano: gian- Gamma-Ray Bursts are bright flashes of high energy photons 
carlo.cusumano@ifc.inaf.it that can last from about 10 milliseconds to 10 minutes. Their 
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origin and nature puzzled the scientific community for about 25 
years until 1997, when the first X-ray afterg low of long (> 2 s 
duration) burs ts were detected ( Costa et al. 1997 ), and the firs t 
optical ( Ivan Paradiis et al. 1997b and radio dFrail et al. 19971) 
counterparts were found. These measurements established that 
long GRBs are typically at high redshift ( z~ 1.6) and are in 
sub-luminous star-forming host galaxies (iBloom et al. 20021) . 
They are likely produced in core-collapse explosions of a 
class of massive stars that give rise to highly relativistic jets 
(the collapsar model; (IMacFadven et al. 200 lb ). Internal inho- 
mogeneities in the velocity field of the relativistic expanding 
flow lead to collisions between fast moving and slow mov- 
ing fluid shells and to t he formation of internal shock waves 
( Ree s &Meszaros 19941) . These shocks are believed to produce 
the observed prompt emission in the form of irregularly shaped 
and spaced pulses of gamma-rays, each pulse corresponding to 
a distinct internal collision. The expansion of the jet outward 
into the circumburst medium is believed to give rise to "ex- 
ternal" shocks, responsible for producing the smoothly fading 
afterglow emission seen in the X-ray, optical and radio bands 
dMeszaros & Rees 1997b. 



The Swift dGehrels et al. 2004 X-ray Telescope (XRT; 
( Burrows et al. 2005bl) ) is providing a growing number of 
unprecedented observations of the early stages of GRB 
afterglow in the 0.2-10 keV X-ray band. The XRT rapid 
(< 2 min) response to the Swift Burst Alert Telescope 
(BAT; dBarfhelmv et al. 2005 ah ) triggers has already led 
to the discovery of rapid early X-ray declines followed 
by the smoother "standard" X-r ay afterglow compo 



( Taglia ferri et al. 2005at Cusumano et al. 2006a ; 



nents 

Barthelmv et al. 20005b!; 
matic flaring in the 



Vaughan et al. 2006al) . dra- 
!-ray light curves of short 



dFox et al. 20051 iBarfhelmv et al. 2005c : 


Campana et al. 2006 ; 


Vaughan et al. 2006b|) and long bursts ( 


Burrows et al. 2005a; 


Romano et al. 2006; Falcone et al. 2006; 


Burrows et al. 2005c: 


Pagani et al. 2006; Morris et al. 2006) and simultaneous peaks 



at the end of the BAT observation (15-350 keV) and at the 
beginning of the XRT observation (0.2-10 keV) of some 
long bursts (GRB 050730, GRB 050820a, GRB 050822; see 
O'Brien et al. 2006). Thanks to its fast response, and precise 
source localization, about 3' in BAT, 3"-5" in XRT and 0.3" 
in UVOT (Roming et al. 2006), Swift is able to quickly alert 
ground-based telescopes to locate the optical counterpart and 
get redshift measurements before the object becomes too faint. 
GRB 050904 triggered the BAT on 2005 September 4 



at 01:51:44 UT dCummings et al. 2005b . The burst was lo- 
cated on-board at RA J2 ooo=OOh54m41s, Decj2ooo=+14° 08' 
17" with an uncertainty of 3' radius (90% confidence level) 
and the spacecraft was quickly pointed towards it. The XRT 
on-board centroiding algorithm failed to detect the counter- 
part due to the presence of a hot CCD detector column. The 
burst was long and bright with duration Tgo- (225+10) s 
and a 15-150 keV fl uence of (5.4+0.2) x 10~ 6 erg crrT 2 
(Sakamot o et al. 20051) . UVOT did not detect the burst coun- 



terpart down t o a 3 sigma upper li mit of about 21 mag in all 
its six filters dCucchiara et al. 2005 ). The optical detection was 



optical flar e was detected during the prompt high-energy emis- 
sion phase (Bo er et al. 2006HGendre et al. 20061) . 

Early ph otometry indicated a high redshift (z>5, 



dReichart 20 05)). A photometric redshift z = 6.1 J. ,L was mea 



sured by Tagliaferri et al. (2005b) and confirmed by a Subar u 
spectroscopic measurement of 6.29+ 0.01 (iKawai et al. 2005 ). 
A break at T/, = 2.6+1.0 days wa s also found in the J-band 
light curve ( Tag liaferri et al. 2005B). 

Such a high redshift means that this explosion happened 
12.8 billion years agcQ, corresponding to a time when the 
Universe was youn g (< lGyr), close to the re-ionization era 
dBeckeretal. 200lh . This gave GRB 050904 the distinction of 
being the most distant cosmic explosi on ever observed: th e 
previous record for a GRB was 4.5 ( Andersen et al. 2000h . 
the most dist ant quasar known is at a redshift of 6.4 
( Fan et al. 20031). and the most distant galaxy is at a redshift 
of ~ 7 dKneib et al. 2004 . 

Here we present the analysis of the BAT and XRT obser- 
vations of GRB 050904. The details on data reduction are in 
Section 2; temporal and spectral analysis results are reported in 
sections 3 and 4, respectively. In section 5 we draw our con- 
clusions. Our results have already been shortly summarized by 
Cusumano et al. (2006b), but this paper contains a complete de- 
scription of the phenomenology of the GRB 050904. Compared 
to the Watson et al. (2005) paper, that already presented the 
XRT data analysis, we add: ( i) the detailed BAT data analysis, 
( ii) complete and correct BAT and XRT light curves, ( Hi) time 
resolved spectral analysis on a finer temporal grid, (iv) spectral 
analysis of BAT and XRT simultaneous data. Hereafter, errors 
are reported for a 90% single parameter confidence level. 

2. Observations and Data reduction 

The BAT data were analyzed using the standard BAT analysis 
software included in the HEAsoft 6.0.4 package and described 
in the Swift BAT Ground Analysis Software Manual! BAT data 
from 306 to 525 s after the burst onset were telemetered in the 
masktag-lc observing mode that accumulates data in only four 
energy bands. As a consequence, no spectral energy analysis 
was performed on BAT data for this time interval. 

GRB 050904 was observed by the XRT from 161 s up to 10 
days after the burst onset, overlapping the BAT observations for 
about 364 s. The XRT observation started before the end of the 
high energy prompt emission. Data were accumulated in WT 
mode up to 573 s after the trigger time, while all the other data 
were obtained in PC mode. In the WT mode only the central 
8' of the field of view is read out, providing one dimensional 
imaging and full spectral capability with a time resolution of 
1.8 ms. The PC mode provides, instead, full spatial and spectral 
resolution with a time resolution of 2.5 s. 

XRT data were calibrated, filtered and screened using the 
XRTDAS package included in the HEAsoft 6.0.4, as described 
in the XRT Software User's Guide 2 . Only observing time inter- 
vals with a CCD temperature below -47 degrees Celsius were 



first made by the robotic observatory TAROT that begun to ob- 
serve GRB050904 only 81 s after the Swift trigger. A bright 



1 We used standard cosmological parameters of O M =0.27, 
fi A =0.73, #o=71 kms-' Mpc" 1 

2 http://swift.gsfc.nasa.gov/docs/swift/analysis/ 
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Rest frame Time (s) 

Fig.l. The 0.2-10 keV light curve of GRB 050904 as observed by the BAT and XRT. The observed BAT count rates were 
extrapolated into the XRT 0.2-10 keV band using a conversion factor evaluated from the BAT best fit spectral model (TableQ~|l. 
The observed XRT count rates were converted into flux Fo.2-10 keV (scale on the right side) using the best fit spectral parameters 
listed in TableQ] The equivalent isotropic luminosity (scale on the left side) was calculated as Fo.2-10 keV 47r d\ (l+z) (r ~ 2) where 
cLl is the luminosity distance and the last term is the K correction factor for emission modeled by a power law. The error bars are 
given by the quadrature sum of the count rate statistical error and the estimated uncertainties in the conversion factors. The inset 
shows BAT and XRT fo.2-10 keV for the first orbit and the best fit model for XRT data. Note that small differences in the the burst 
history shown in the inset (Fo.2-10 keV units) and the main picture, which shows the evolution in luminosity units, are due to the 
photon index dependence in the K correction factor and its measurement in discrete time intervals. 



used. The total exposures after all the cleaning procedures were 
2.8 ks and 127.3 ks for data accumulated in the WT and PC 
modes, respectively. We used a 0-2 and 0-12 grade selection 
for data in the WT and PC modes, respectively. Such a selec- 
tion provides the best combination of spectral resolution and 
detection efficiency. The GRB was imaged far from the CCD 
hot columns and no corrections for hot columns inside the pho- 
ton extraction regions was necessary. 

T he position of the burst was refined by on-ground anal- 
ysis ([Palmer et al. 2 005). The BAT burst position is RAj2ooo = 
00h54m53s, Dec J2 ooo = +14°04'52 ", with an uncertainty of 2!6. 
This is 3>'.9 from the o n-board position and 0f5 4 from the near- 
IR afterglow position dNvsewanderet al. 20051) . The XRT after- 
glow position derived with xr tcentroid (vO.2.7) a nd including 
the latest boresight correction (IMoretti et al. 20061) is RAj2ooo = 
00h54m50s.8, Dec J20 oo = +14°05'09'.'0, with an uncertainty of 
3'.'5. The XRT derived coordinates are 35'.'9 from the BAT ones 



and 074 from the near-IR counterpart dNvsewanderet al. 20051) . 

The BAT and XRT times are referred to the GRB 050904 
onset T =2005 Sep 4, 01:51:44.3 UT. 



For the measured redshift z = 6.29, the 15-350 keV BAT 
band corresponds to a 109-2551 keV band in the burst rest 
frame while the 0.2-10 keV XRT band corresponds to a 1.4- 
73 keV band. The observed timescales are stretched by a factor 
(1 + z) with respect to the rest frame ones. In the following 
the GRB phenomenology is presented and discussed from the 
point of view of the source rest frame. 

3. Timing Analysis 

The timing analysis for the XRT data was performed by select- 
ing events from a region centered on the source with a radius 
of 20 (47.2) and 35 (82.6) pixels (arcsec) for WT and PC data, 
respectively. The background was estimated from regions suf- 
ficiently offset (> 2 arcmin) from the source position to avoid 
contamination from the PSF wings and free from contamina- 
tion by other sources. 

The intensity of the source caused pile-up in the PC data 
up to 8 ks from the burst onset. The pile-up correction was 
performed by excluding a central region of 4 pixels radius and 
dividing the extracted count rate by the fraction of lost point 
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spread function (PSF; 58%). For the rest of the observation we 
used the full circular extraction region. 

WT data were also corrected for the fraction of PSF not 
included in the extraction region (7%). 

FigureQ]shows the evolution of the GRB flux and luminos- 
ity. The BAT light curve displays three main peaks: two short 
peaks (~2 s long) at T+3.8 and at T+9.3 s, and a main long- 
lasting peak at ~T+13.7 s, where T is the time of the burst 
onset. Emission in the BAT energy range continues up to al- 
most T+77.7 s with a weak peak at ~T+65 s, coincident with 
the first peak of the XRT light curve. The BAT and XRT light 
curves overlap between T+23 and T+69 s. The early XRT light 
curve shows a steep decay with a slope a = -2.07 + 0.03 with 
three flares superimposed at T+65 s, T+126 s and T+171 s. 
These flares can be modeled by a linear rise lasting 26.6, 5.3 
and 4.7 s, plus an exponential decay with decay time of 4.5, 
10.98 and 5.2 s, respectively. The best fit model of the first 
orbit of the XRT light curve is shown in the inset of Fig. [T] 
Although interrupted by observing constraints imposed by the 
Swift orbit, the light curve from GRB 050904 reveals highly ir- 
regular rate variations likely due to the presence of flares up to 
T+1.5 hours. At later times the flaring activity is not detected 
and only a residual emission, 10 5 times lower than the initial 
intensity, is visible. Note that the XRT light curve presented in 
Watson et al. (2005) shows an evident discontinuity at the end 
of the WT observation segment likely caused by an error in his 
flux conversions for the following PC mode data, all of which 
are systematically too low by about a factor of four. 

Figure [2] shows the rest frame XRT light curves of 
GRB 050904 in the 1.4-14 keV and 14-73 keV bands (top and 
middle panels). The two light curves are binned before conver- 
sion to the rest frame in order to have at least 40 counts per 
bin in both bands. The hardness ratio (bottom panel) is defined 
as H/S, where H (hard) and S (soft) are the high energy and 
the low energy bands. The H/S plot reveals a significant shift 
to softer energies with time during the first 80 s, with the ex- 
ception of the flaring episode at T+65 s, where the hardness 
ratio peaks, too. Later emission shows no evidence of further 
softening until the tail of the last flare around T+5500 s. 

4. Spectral Analysis 

The spectral analysis of GRB 050904 was performed by se- 
lecting two sets of time intervals for the BAT and XRT ob- 
servations, corresponding to characteristic phases of the light 
curve evolution. The BAT spectra were accumulated in the 14- 
150 keV observed band in six time intervals up to 41.6 s from 
the burst onset (Table 1). No BAT spectrum was accumulated 
from T+41 .6 to T+72 s because BAT data were in the masktag- 
lc observing mode with only four energy bands. The XRT spec- 
tra were accumulated in thirteen time intervals from T+23 s 
to T+8173 s from the burst onset (Table 1). The WT and PC 
spectra were extracted from the same regions used for the tim- 
ing analysis. Instrumental energy channels below 0.3 keV and 
above 10 keV for PC and WT spectra were ignored and the 
background was evaluated in regions sufficiently offset (> 2 ar- 
cmin) from the source position to avoid contamination from the 
PSF wings and free of contamination from other sources in the 




100 _ 1000 
Rest frame Time (s) 



Fig. 2. Hardness ratio evolution of GRB 050904. The top and 
middle panels show the count rate evolution of XRT data in 
two different energy bands in the rest frame. 



field of view. Moreover, energy channels between 0.5 and 0.6 
keV were excluded because of the presence of a sistematics in 
such a range due to a time dependent energy scale problem at 
low energjlj a sistematics in the gain offset produces a negative 
residual in the 0.5 keV region increasing in the fit procedure the 
N H estimation. The BAT spectra were modeled with a power 
law with photon index T (F(E) oc £ r+1 ) while the XRT spectra 
were modeled with a power law plus two absorption compo- 
nents: one for the intrinsic absorption in the host galaxy and 
one for the Galactic absorption. The latter was fixed to the line - 
of-sight value of 4.93xl0 20 cirr 2 jDickev & Lockman 199dh . 
The model gave acceptable x\ values for all the selected inter- 
vals (Table 1). Mo re complex models, such as a Band function 
(IBand et al. 1993b . cannot be constrained by the data. Figure[3] 
shows the evolution with time of the photon index T. The BAT 
spectra have F ~ -1.2, consistent w ith typical values of the 
ffBand parameter of the Band model JPreece et al. 20001) . This 
strongly suggests that the BAT observes the low energy part of 
the Band function and that the peak energy of the GRB spec- 
trum is above 150 x (1 + z) keV in the source rest frame. If we 
exclude the spectrum of the first XRT flare at T+65 s, the XRT 
photon indices show a clear decreasing trend from about -1.2 
to about -1.9 in the first T+200 s. No further spectral evolution 
is present in later XRT data, in agreement with the hardness ra- 
tio curve. The BAT and XRT photon indices are in good agree- 
ment in the overlapping region. We also detected in the XRT 
WT spectra (T < T+67.1 s), a highly significant absorption in 
excess of the Galactic value. The intrinsic absorption column 
decreases with time at high significant level. In the time inter- 
vals where the burst was observed in PC mode (T > T+67.1 
s), due to the decreased statistics, only upper limits were mea- 
sured. Table Q] shows the best fit parameters for each of the 
selected time intervals. 

The GRB 050904 spectral evolution is also quite evident 
in Fig. [4] where the £ 2 F(E) deconvolved spectra (equivalent to 
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Interval 


Time (s) 


N H (10 22 cm 2 ) 


r 


xi (dof) 


Flux (10 1J erg cm 2 s ') 




Start 


Stop 






0.2-10 keV 


15-350 keV 




BAT 1 


-1.43 


2.69 


— 


-1.2 +0.4 


1.2 (57) 


1.4 


22.9 


2 


2.69 


4.89 


— 


-1.05±0.16 


0.86 (57) 


3.2 


90.8 


3 


4.89 


10.1 


— 


-1.36±0.21 


0.97 (57) 


3.4 


30.9 


4 


10.1 


20.4 


— 


-1.17±0.08 


0.95 (57) 


3.6 


66.8 


5 


20.4 


30.6 


— 


-1.22±0.10 


0.93 (57) 


3.0 


45.7 


6 


30.6 


41.6 


— 


-1.5 ±0.3 


0.88 (57) 


2.0 


9.9 


XRT 1 


23.2 


28.7 


5.73±4.2 


-1.19±0.1 


0.77 (62) 


3.5 


— 


2 


28.7 


36.9 


5.5±2.5 


-1.34±0.08 


0.98 (95) 


2.5 


- 


3 


36.9 


50.6 


3.4±2.2 


-1.33±0.08 


0.78 (89) 


1.3 


— 


4 


50.6 


58.8 


7.7 ±4.5 


-1.85±0.1 


1.12 (56) 


1.4 




5 


58.8 


67.1 


4.2±2.0 


i ct\ > f\ r\c\ 

-1.50±0.09 


1.14 (73) 


1.7 




6 


67.1 


79.8 


1.5+1.4 


-1.86+0.13 


0.94 (37) 


0.54 




7 


79.8 


159.4 


<6.4 


-1.80+0.15 


1.12 (23) 


0.12 




8 


159.4 


244.4 


<6.4 


-1.97±0.24 


0.90 (7) 


0.05 




9 


628 


848 


<5.2 


-1.80±0.24 


0.92 (7) 


0.02 




10 


848 


1040 


<6.8 


-1.86±0.14 


0.90 (35) 


0.08 




11 


1452 


1863 


<5.8 


-2.01 ±0.22 


0.80 (17) 


0.02 




12 


2275 


2618 


<6.9 


-1.90±0.14 


1.26 (47) 


0.04 




13 


3045 


8173 


<4.0 


-1.97±0.12 


1.24 (35) 


0.008 





Table 1. BAT and XRT spectral analysis results. The BAT fluxes in the XRT band are extrapolated from the best fit models. 
Quoted errors are at the 90% confidence level. 



vF(v) spectra) selected in 5 ad hoc intervals (bottom panel in 
Fig. 3) are shown together with the best fit spectral models. The 
first spectrum (A) is accumulated in the time interval T-1.43 
s - T+23.2 s, when only BAT observed the burst. The second 
spectrum (B) is from T+23.2 s to T+41.7 s when the GRB is 
simultaneously observed by BAT and the XRT. The following 
spectra (C, D, E) refer to the time intervals T+41.6 s - T+82.3 
s, T+82.3 s - T+224 s and T+628 s - T+8163, respectively, 
when only the XRT detected the burst. The energy distribution 
is clearly softening with time starting from interval C. The fit in 
the broad energy band of interval B is also well modeled by a n 
absorbed power law while a Band function dBand et al. 1993b . 
cannot be constrained by the data. Table [2] shows the best fit 
parameters for each of the selected time intervals. 

We also evaluated the contribution to the total fluence in the 
1.4-73keV band of the three flares (T+65, T+126 and T+171 
s) superimposed on the early XRT light curve. The fluence val- 
ues over the continuum are ( 1 .2 + 0.08) x 1 0~ 6 , (4.7 + 0.5) x 1 -8 
and (5.8 + 0.6) x 10~ 8 ergcrrT 2 , respectively. The fluence of the 
XRT continuum over the first orbit (i.e. from 23.2 to 244.4 s) is 
(4.9 + 0.3) x 10~ 6 ergcirr 2 . The extrapolated 1.4-73keV BAT 
fluence in the time interval from the burst onset to the start 
of the XRT observation is (4.1 + 0.2) x lO^ergctrT 2 . The 
three XRT flares are therefore 5%, 1% and 1% of the total 
1.4-73 keV emission observed up to T+244 s. The 1.4-73 keV 
fluence in the remaining part of the XRT observation is 1.8 x 
10 -6 ergctrT 2 . This value is only a lower limit because of the 
observing gaps. 



5. Discussion 

GRB 050904 was a long, multi-peaked, bright GRB 
with strong X-ray flaring activity lasting up to 1- 



-0.8 
-1.0 
S -1.2 

I -1.4 

o 

-1.6 

D. 

-1.8 
-2.0 
-2.2 
10 5 ' 

1 10" 

1 1 ° 4 ° 

o 
_- 

I 10" 

10" 



BAT 

XRT-WT 
XRT-PC 




10 100 

Rest frame Time (s) 



• t 

E 



Fig. 3. Spectral evolution of GRB 050904. The top panel illus- 
trates how the photon index T of GRB 050904 changes during 
the observation. In the bottom panel the burst evolution is plot- 
ted to show how the time intervals for spectral analysis were 
selected. Vertical bars indicate the time intervals selected for 
BAT and XRT spectral analysis whose results are reported in 
Table 2 and Fig. 4. 



2 hours in the source rest frame (or 5 — 6 xlO s 
in the observer frame). X-ray flares are detected in 
nearly half of the S wift long bursts (e .g. GRB 050406 
Burrows et al. 2005at [Romano et al. 2006I) . GRB 050502 
Falcone et al. 



(Burrows et al. 2005at iFalcone et al. 2006I). GRB 50607 
dPagani et al. 2006I) . GRB 0507 13A dMorris et al. 2006k The 
rise time and decay time of these flares seen at lower z are 
frequently very fast with a ratio between the duration and peak 
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Interval 


Time (s) 


N H (10 22 cm.- 2 ) 


r 


Xi (dof) 


Flux (10-" 


erg cm 2 s 1 ) 






Start 


Stop 








0.2-10 keV 


15-350 keV 


BAT 


A 


-1.43 


23.2 




-1.19 ±0.07 


0.96 (57) 


3.1 


52.8 


BAT+XRT 


B 


23.2 


41.6 


7.5+1.5 


-1.37±0.02 


1.0(189) 


2.6 


11.9 


XRT 


C 


41.6 


82.3 


4.7+1.2 


-1.63±0.05 


1.19 (205) 


1.33 




XRT 


D 


82.3 


224 


< 2.1 


-1.7±0.1 


1.24 (22) 


0.03 




XRT 


E 


628 


8163 


<2.9 


-1.74+0.07 


1.243 (80) 


0.006 





Table 2. BAT and XRT spectral analysis results of relative to 5 ad hoc selected intervals of the GRB evolution (see bottom panel 
in Fig. 3). The BAT fluxes in the XRT band are extrapolated from the best fit models. Quoted errors and upper limits are at the 
90% confidence level. 




energy (keV) 

Fig. 4. Spectral evolution of GRB 050904; BAT and XRT 
£ 2 F(E) spectra (equivalent to vF(v) spectra) deconvolved from 
the detector response for five selected time intervals (see bot- 
tom panel in Fig. 3). 



time 6t/t <K 1. These features make difficult to explain these 
flares with mechanisms associated with the external shock 
(Zhang et al. 2006) 

The variability of GRB 050904 X-ray light curve is even 
more dramatic than the typical Swift afterglow showing flares, 
and the amplitude and rise/fall times of the flares are consis- 
tent with the behavior of nearby (z <1) long GRBs (Fishman 
& Meegan 1995). This suggests to interpret them as late 
internal shocks related to central engine activity. In this 
scenario they would have the same or i gin as the prompt 
gamma-ra y emissio n dZhang et al. 20061: iNousek et al. 2006; 
O'Brien et al. 20061) . This would require that the central en- 
gine remain acti ve up to at least 5000 s, consistent with the 
collapsar model (iMacFadyen et al. 200 ll) . which allows central 
engine activity for up to a few hours. Then, BAT and XRT have 
likely recorded a single continuous observation of long lasting 
prompt emission where the standard X-ray afterglow compo- 
nent is hardly detectable because its contribution is drowned 
by the intense and long lasting activity of the central engine. 

TAROT ob se rved a flare i n the optical band 
(Boe ret al. 20061: Idendre et al. 20061) which is simultane- 
ous with the X-ray flare at T+65 s . Later ( > T+1200 s) 



times in X-rays is indeed present. The multiwavelength fit 
by Tagliaferri et al. (2005b) suggests that the optical data 
taken after T+1200 s are in agreement with standard afterglow 
emission and jet lateral spreading at about T+80 ks. Physical 
interpretations to explain the simultaneous X-ray and optical 
flare at T+65 s are extensively discussed in Gendre et al. 
(2006). They showed that delayed external shock from a thick 
shell (Piro et al. 2005), synchrotron emission from reverse 
shock (Fan & Wei 2005), inverse Compton emission from 
reverse shock (Kobayashi et al. 2005) and inverse Compton 
from late internal shock (Wei et al. 2006) cannot satisfactory 
explain the simultaneous X-ray and optical flare. On the other 
hand, the late internal shock model proposed by Zou et al. 
(2006) could account for the broad band data of the first flare. 

In the time interval from T+23 to T + 244 s, the ob- 
served intensity underlying the XRT flares decays as f 
with a = -2.07+0.03. An initial steep decay of the X- 
ray emission has been observed in many other GRBs de- 
tected by Swift dTagliaferri et al. 2005at iNouseket al. 2006 : 
O'Brien et al. 20061) . The decay slope together with the XRT 



optical and infrared observations (Tag liaferri et al. 20 05b) 
are too sparse to test if the flaring activity observed at later 



energy index j3 — T + 1 ~ -0.2 measured up T+50 s are in good 
agreement with the interpretatio n of the observed emissio n as 
due to high-latitude emission (iKumar & Panaitescu 20001) . In 
this model, the tail of a peak is the emission from the shocked 
gas moving at angle 6 > 1/y (where y is the Lorentz factor) 
relative to the observer. The higher the angle 6, the later the 
photon arrival-time and the weaker the relativistic beaming of 
the emission, leading to a f decay with a — (3 - 2. After T+50 
s, due to the decrease of (3 to about -0.7, the predicted slope 
would be steeper than the measured -2. This deviation could 
be reconciled with the high-latitude emission assuming that the 
delayed radiation from the outer parts of the emitting curved 
shell is softer and b righter than the radiation along our line of 
sight dKumar et al. 20061) . 

The decrease of the photon index around T+50 s could be 
interpreted as an indication of a shift of the spectral peak en- 
ergy (E p ) towards lower energies, but poor statistics and the 
narrowness of the XRT energy range do not allow us to verify 
this hypothesis. 

We detected highly significant absorption in excess of the 
Galactic value. The host column density decreases with time 
at a high significance level. This is consistent with the idea 
that the circum-burst absorbing material is photoionized by 
the high-energy photons of the burst (Perna & Loeb 1998). 
Evidence for such a decrease has been found for GRB 980329 
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dFrontera et al. 20001) and GRB000528 dFrontera et al. 2000h 
but neither of them has a comparable data quality. We do not 
confirm evidence found by Watson et al. (2006) of an increase 
of the column density at the peak of the first X-ray flare. 

Our lack of knowledge concerning the peak energy of the 
BAT and XRT spectra does not allow a precise estimate of the 
total energy released by GRB 050904. However, we can calcu- 
late lower and upper limits to the isotropic-equivalent radiated 
energy up to 244 s from the burst onset, i.e. including con- 
tributions from the first three XRT flares. To evaluate the lower 
limit to Ei S0 we integrated the best fit power law spectral en- 
ergy distributions in the (1 - 200)x(l+z) keV band and in the 
(1 - 10)x(l+z) keV band for the BAT and XRT, respectively. 
The standard energy range 1 — 10 4 keV (rest frame) was used 
to evaluate the upper limit to £j so . We obtained 6.6 x 10 53 erg 
< Eiso < 3.2 X 10 54 erg. Additional contributions from the later 
flare portions are only a few percent in both the upper and lower 
limit. The large isotropic equivale nt energy of this b urst is in 
agreement with the Amati relation ( Amati et al. 20021) with an 



E p of about 1500 keV in the rest frame. This is consistent with 
our non-detection of a peak energy in BAT spectral fit. 

The break observed in the optical and infrared afterglows 
( Tagl iaferri et al. 20 05b) at Tb = 2. 6± 1.0 days (observer frame) 
and the range of E mo evaluated above imply a jet half-opening 
angle §j et between 2° and 4°, assuming a radiative efficiency 
r] = 20% and a circumburst medium density n = 3 cirT 3 . 
This angle estimate is consistent with those obtained by mod- 
elling optical light-curve breaks observed in 10 pie-Swift GRBs 
(Panaitescu 2005). It corresponds to a collimation-corrected 
energy E y between 0.4 and4.4xl0 51 erg. This is well within the 



E y distribution of GRBs with known redshift (Fra il et al. 2001 



Bloom et al. 2003). Consistency with the Ghirlanda relation 



( Ghirlan da et al. 20041) constrains the rest frame peak energy 
of the average spectrum to be between 560 and 1300 keV. 

In the optical afterglow, the measured pre-break slope de- 
cay a a = -0.7 + 0.2 (Tagliaferri et al. 2005) requires an electron 
index of p = (4/3)a c +l ~ 1.9±0. 3 in a local uniform interstellar 
medium. The expected intrinsic optical spectral slope should 

be p o = (p-l)/2 0.45 ±0.15 while the observed one is fi = 

- 1 .25+0. 15 (Tagliaferri et al. 2005). This discrepancy could be 
naturally explained assuming the presence of a bit of dust (Ay = 
( 1 .25-0.45)/( 1 +z) = 0. 1 ) in the local host medium. In the X-ray 
regime the slowest decay of a possible underlying continuum, 
inferred by fitting the lowest points in Fig. 1 after T+1000 s, 
is a x ~ -1.2. This is marginally consistent with the slope ex- 
pected for the forward shock emission for the value of electron 
index inferred from optical data: a x = (3p-2)/4 ~ 0.95+0.20, 
under the hypothesis that the cooling frequency is between op- 
tical and X-rays. The expected spectral index in X-rays is f$ x = 
p/2 — 1.0+0.1, consistent with the observed/? v •= - 1.0+0.2 af- 
ter T+ 1000s. The continuum component inferred by the fitting 
the the lowest points in Fig. 1 is, therefore, marginally consis- 
tent with a standard forward shock emission and the properties 
of the optical emission reported by Tagliaferri et al. (2005). 
However, we cannot prove that the estimated underlying X-ray 
continuum and the late optical light curve are the same forward 
shock emission because the observational gaps do not allow to 
see the minima of the X-ray light curve. 



Figure 5 shows how GRB 050904 would appear in the Swift 
0.2-10 keV band if it were at redshifts other than z=6.29. The 
observed intensity varies inversely with the square of the (1+z) 
factor. A further dependence on z is due to the K correction 
that accounts for the redshift dependence of the luminosity in a 
given wavelength band. The observed timescale undergoes dif- 
ferent stretching factors with respect to the rest frame: the burst 
would appear longer at higher redshift. The horizontal line in 
Fig. 5 gives an indication of the sensitivity limit of BAT, while 
the vertical line marks the start of the XRT follow-up. Starting 
from a redshift lower than ~ 2, BAT would have observed all 
later flares. GRB 050904 would have been detected up to a dis- 
tance corresponding to z~10. Comparing the intensity of other 
GRBs observed by Swift (see Fig. 3 in O'Brien et al. 2006) with 
the flux that we would observe if GRB 050904 had exploded 
at a redshift lower than 2 ( ~ the average redshift of the Swift 
GRBs) we see that this burst was intrinsically bright in X-rays. 

After one year of operations, Swift has detected four con- 
firmed high-redshift ( it 4) GRBs (out of 24 with known red- 
shifts). Figure 6 shows the K-corrected 0.2-10 keV luminos- 
ity evolution of these GRBs, including both BAT and XRT 
data set. BAT light curves are obtained by extrapolation of the 
15 - 150 keV light curve to the XRT energy band. In addition 
to GRB 050904 this sample includes GRB 050730 at z=3.969 
(Chen et al. 2005), GRB 050505 at z=4.27 (Berger et al. 2006) 
and GRB 050814 at z=5.3 (Jakobsson et al. 2006). All of the 
bursts are exceptionally luminous and long-lasting (as mea- 
sured in the source rest-frame) and are among the brightest 
GRBs ever observed. Their exceptional intensity is not due to 
selection effects, since the 15 - 150 keV burst fluxes are well 
above the BAT detection threshold. Rather, their unique prop- 
erties could be due to their likely origin within some of the first 
star-forming regions in the Universe; Woosley & Heger (2006) 
suggested that the very low metallicities of the progenitors at 
these epochs may provide an explanation. A more reliable con- 
clusion about systematic differences and similarities in lumi- 
nosities and durations of high redshift GRBs will require an 
increase in the sample size, which should come in future years 
of Swift operation. 

Detecting high-redshift GRBs with Swift, and measuring 
their redshifts with ground-based spectroscopy, is of substantial 
interest because of the link between long-duration GRBs and 
the star formation rate. The GRBs with measured redshift can 
be used to infer the cosmological star formation history, with 
relatively minor (or in any case unique) selection effects by 
comparison to other methods (Porciani & Madau 2001; Lamb 
& Reichart 2000; Bromm & Loeb 2002; Natarajan et al. 2005). 
A preliminary estimate of the star formation rate derived from 
Swift bursts (Price et al. 2006) shows, within current broad un- 
certainties, a flat or (at the highest redshifts) slowly-declining 
star formation rate, consistent with results obtained from color- 
selected galaxy observations (Bunker et al. 2004). 

6. Summary and conclusion 

We have presented the results of the analysis of BAT and XRT 
observations of the high redshift (z=6.29) GRB 050904. The 
GRB light curve is characterized by the presence of a long flar- 
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Fig. 5. GRB 050904 light curves in the 0.2 - 10 keV band as 
they would be observed if the burst were located at different 
redshifts. Each light curve includes both BAT and XRT extrap- 
olations. The horizontal line gives an indication of the sensitiv- 
ity limit of BAT, while the vertical line marks the start of the 
XRT follow-up in the observing frame. 
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Fig. 6. The K-corrected 0.2-10 keV time histories of the high- 
redshift ( i 4) GRBs observed by Swift. Each light curve in- 
cludes both XRT and BAT data extrapolated to the XRT energy 
range. The K correction was performed for an average photon 
index for simplicity. 



ing activity lasting up to 1-2 hours after the burst onset in the 
burst rest frame, with no evidence of a smooth power-law decay 
following the prompt emission as seen in other GRBs. We inter- 
pret the overall phenomenology of GRB 050904 as long lasting 
prompt emission where the standard X-ray afterglow compo- 
nent is hardly detectable because its contribution is drowned by 
the intense and long lasting activity of the central engine. The 
spectral energy distribution softens with time, with the pho- 
ton index decreasing from -1.2 during the BAT observation to 
-1.9 at the end of the XRT observation. The dips of the late X- 
ray flares may be consistent with an underlying X-ray emission 
arising from the forward shock and with the properties of the 
optical afterglow reported by Tagliaferri et al. (2005b). 



Highly significant absorption in excess of the Galactic 
value has been detected. The intrinsic hydrogen-equivalent col- 
umn density shows a significant decreasing with time that we 
interpret as due to the photoionization of the circum-burst ab- 
sorbing material by the high-energy photons of the burst caus- 
ing a graduale reduction of the opacity. 

We have calculated lower and upper limits to the isotropic - 
equivalent radiated energy Ej SO up to 244 s from the burst onset, 
i.e. including contributions from the first three XRT flares. We 
obtained 6.6 x 10 53 erg < E iso < 3.2 x 10 54 erg. This range 
of E, w and the break observed in the optical and infrared af- 
terglows imply a jet half-opening angle i9] et between 2° and 
4°, assuming a radiative efficiency rj = 20% and a circumburst 
medium density n = 3 cirT 3 . 
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